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Abstract. We establish a variational framework for nonlinear instabilities in a setting of the 
ideal magnetohydrodynamic (MHD) equations. We apply a variational method to various kind 
of smooth steady states which are shown to be nonlinearly unstable for both incompressible 
and compressible ideal MHD equations. Destabilizing effect of compressibility is justified as 
well as stabilizing effect of magnetic field lines arising in MHD dynamics, which distinguishes 
from the Rayleigh- Taylor instability in the absence of magnetic field lines. 



1. Introduction 

Magnetohydrodynamic equations (MHD) serve as an important model for fluid and gas 
dynamics and hydromagnetic instability is a fundamental phenomenon in nature, for instance, 
oceans, atmosphere, and plasma. MHD instability possesses extensive applications in both 
laboratory plasmas and astrophysics such as in nuclear fusions, compression of thin foils for X- 
ray production, and stellar dynamics. Nevertheless, there have been not many analytical results 
to date due to its structural complexity such as the presence of shock waves. Furthermore, 
nothing has been known about nonlinear instabilities for MHD equations despite its importance 
and variety of instabilities. 

The main purpose of this article is to present a variational framework in the passage from 
linear to nonlinear instability in a setting of the ideal MHD system and derive nonlinear in- 
stability around different steady states for both incompressible and compressible ideal MHD 
equations. 

We consider the equations of ideal magnetohydrodynamics(MHD) for inviscid flows: 



Pt + V-(pV) = 0, 
p{V t + V -VV) = {V xB)xB-VP + pg ) (1.1) 
B t = Vx(V x B), 
V-B = 0. 

Here p is plasma density, V velocity, B magnetic field, P plasma pressure, and g is the 
gravitational field. The x axis is taken along the gravitational field, which is assumed to be 
uniform: 



9 = (9,0,0) 
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4:(Po + Wo)=Po9- (1-2) 



The condition for steady state (p , S , P ) with Vq = is then 
We assume for density profile that 

min p Q > c > 0. (1.3) 

X 

Our domain is t > and 

D = {0 < x < 2vr, < y < 2tt} . 
We assume periodic conditions at the boundary for V and 5. Both incompressible and com- 
pressible fluids are considered: In the incompressible case, 

V-V = 0. 

In the compressible case, 7 is an adiabatic index which relates the pressure P to the density p 
by 

P(p) = Cp\ (1.4) 

— = C 7 p^ 2 Vp =: q (p) Vp. 
P 

Or equivalently 

p t + v ■ Vp + 7PoV • w = (1.5) 

We now consider perturbations (<x, i>, S, p) around such a steady state (p , So, Vo = 0, So) of the 
form: 

Vl = Vl (t, x, y) , w 2 = ^2 (t, x, y) , v 3 = v 3 (t, x, y) , 
Si = B x (t, x, y), B 2 = S 2 (t, x, y) , S 3 = S 3 (t, x, y) . 
Equations for perturbed quantities take the form: 

°t + V • [(p + a) v] = 0, 
(p +a)(v t + v- Vv) = (V x B ) x S + (V x S) x S (1.6) 

+ (Vx5)xB-Vj) + 4 
S t = V x (v x (S + 5)) , 
V- B = 0. 
We also obtain its linearized system: 

<rt + V-(p v) = 0, (1.7) 
p v t = (V x S ) x S + (V x S) x S - Vp + of, 
£ t = V x (u x S ) , 
V-S = 0. 

In the incompressible case, v also satisfies, in both linear and nonlinear system, 
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It is crucial to make an equivalent second-order linearized system in one quantity v and use it 
through our variational formulation: By taking t-derivative of f-equation in (jl.7|) and plugging 
a- and -B-equation in (jl.7J) (also ()1.5|) in the compressible case) into the resulting equation, we 
obtain the following second-order linear hyperbolic PDE for velocity v: 



p v tt = (Q-V)B + (B -V)Q-V(B -Q+ Pt )-gV ■ (p v) =: L (v) , (1.8) 

where Q = V x (v x B ) and we have used (jl.llj) . 
For notational convenience, 

Notation 1. For any u and v G L 2 (0, 2-71") , 

< u,v >= / u-vdx, (u,v) = / p u-vdx, 
Jo Jo 

For any u and v G L 2 (D) , 

< u, v >= J J u ■ v dxdy, (u, v) = J J p u ■ v dxdy. 

Notation 2. ||/|| = (f,/) 1 ^ 2 , \\f\\ s = {J2 a ll^ a /l| 2 ) ^ > where a is a multi-index with \a\ < s. 

We make the following variational formulation and obtain a discrete set of eigenvalues for 
the linear operator (J1.8)) : for each wave number k G N, 

\2 <L(u),u > 

X k = sup r , (1.9) 

ui6H K (0,27r), 
u 2 &H K - 1 {0,2rc) 

where d y is replaced with multiplying by k in L (u) and k = 1 or 2. Notice that the RHS of 
(jl.9|) is indeed a function of k. 

In the incompressible case, 112 is replaced by —\u\ x from divergence-free condition for the 
normal modes and (jl.9p reduces to a variational problem for u\ alone. On the other hand, the 
compressible case may not be simplified to a formula for u\ alone. 

Key step of this article is to show that this discrete set {Afc} fegN of eigenvalues characterizes 
the continuum spectral radius of the whole linearized operator by taking the limit as k —>■ 00 
to obtain the bounded least upper bound A > 0: Let 

o o ffr> H (v , v x ) dxdy 
\im\ 2 k = A 2 = sup ^ — -> (1-10) 

k v£H«-i(D) [V,V) 

where coefficients of H depend on p ,Po, B . Then we obtain the following theorem: 
Theorem 1. Let (cr,v) be a solution to \1. 7\ ) and let A 2 > ; then we have 

\\a (t) ,v(t),B (t) I < Ce At \\a (0) , v (0) , B (0) || s+2 , 
where C = C (A, p Q , s) . 
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Thanks to this theorem, we can locate a dominant eigenvalue and pass to nonlinear instability. 
We establish the following dynamical instability for the fully nonlinear MHD system around 
different steady states satisfying A 2 > as in Section 4: 

Theorem 2. Steady states ^p , 0, Bq ~\ of hl.l\) in M. 6 <3jl satisfying A 2 > are indeed nonlinearly 

unstable: For any s large, there exists Eq > 0, such that for any small 5 > 0, there exists a 
family of classical solutions (p s (t, x, y), V 5 (t, x, y), B s (t, x, y)) to (ji.i)) such that 

||/(0, •) - P (OHfl.p) + 11^(0, -)\\ HS{D) + \\B 5 (0, •) - B (•) \\ HS(D) < S, 
but for T s = 0(\\n8\), 

sup s {\\p 5 (t, •) - p (-)|| L1(D) + \\V\t, -)\\ LHD) + \\B S (t, ■) - B (0Li (D) } > eo. 

Remark 1. This indicates the Kruskal-Schwarzs child instability for incompressible ideal MHD 
fluids as 7 — > oo formally. 

Remark 2. The instability time T 5 occurs before the possible blow-up time which is shown in 
the proof and we measure instability in L 1 . 

Rayleigh- Taylor instability is well known as gravity-driven instability in fluids when heavy 
fluid is on top of light one. Linear instability for an incompressible fluid was first introduced 
by Rayleigh in 1883 JB]- Study on linear instability for incompressible ideal MHD system 
in the presence of magnetic field lines could be found in jTH], 0] for instance, which includes 
the classical Kruskal-Schwarzschild instability in the presence of a magnetic field orthogonal 
to the gravitational force. In this case, instability criterion ()4.4j) and the growth rates ([4.2)1 . 
()4.3|) turn out to be the same as the ones in the Rayleigh- Taylor instability as in jTU] without 
effect of magnetic field lines on the instability. On the other hand, when we consider effect of 
a magnetic field parallel to direction of the gravitational force, stabilizing effect of magnetic 
field lines appears as in 1)4.13)1 and 1)4.14)1 . In the presence of a vertical magnetic field, the 
rigidity produced by the magnetic field lines hinders its own way to instability and makes the 
growth rate slower. Condition for linear instability for a compressible fluid in the absence of 
a magnetic field was first derived by Schwarzschild in and since been discussed by many 
other physicists for a certain class of steady states jlj , , [7] , [H] . The full consideration of 
gravity, magnetic field lines and compressibility has been also largely discussed for its linear 
instability by many physicists such as in PQ,|2],5Z1 5 which exhibit very interesting phenomena 
although it accelerates analytical difficulties. 

Despite extensive research and interest in this subject from physical point of view, little has 
been done from mathematical perspective for the MHD system. In addition, the passage from 
linear to nonlinear instability in a conservative PDE system is quite difficult because of the 
following two main obstacles: (1) presence of the continuum linear spectrum and (2) severely 
unbounded high-order terms in PDE systems. No systematic framework has been built up for 
this problem although there have been works towards this subject for specific physical systems, 
for example [H| , [H| , JO] . Variational approach was first introduced by Guo and Hwang JUj i n 
the case of dynamical Rayleigh- Taylor instability for incompressible Euler fluids. However it 
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is not obvious whether it can be extended to MHD instabilities for compressible as well as 
incompressible fluids since MHD has more complicated structure in addition to the analytical 
difficulties coming from compressibility. 

Crucial point is whether and how to locate a dominant eigenvalue in the complex linear 
spectrum of the MHD system. We use the MHD energy principle in order to estimate a sharp 
spectral radius and make extensive use of the variational structure of the linearized MHD sys- 
tem, resulting in more precise and optimal estimates. The radius (|1.1U|) of continuum spectrum 
is obtained as the least upper bound for a discrete set of eigenvalues (|1.9|) of normal growing 
modes and this method is explicit and constructive. We consider three different magnetic fields 
as steady states which give rise to different outcome in the growth rates (|4.3|) . (|4.9|) . (|4.14|) of 
instabilities and different ranges of admissible density profiles for instabilities. We justify stabi- 
lizing and destabilizing effects of magnetic filed lines and compressibility as expected physically. 
Furthermore incompressible case can also be viewed as the limiting case of compressible case 
as 7 — > oo. The article is organized as follows. 

We formulate the variational problems (jl.9|) in Section 2 and show the existence of smooth 
maximizers satisfying the corresponding Sturm-Liouville equations. We then derive a sharp 
growth rate A > for the whole linear system in Section 3 as in the Theorem 1. In Section 
4, we give different examples of steady states which result in different growth rates, Sturm- 
Liouville equations, and different admissibility for density profiles for instability. For instance, 
if a steady magnetic field is parallel to the gravity, we have the fourth-order Sturm-Liouville 
equation ()4.12|) . The compressible case ()4.8|) is more complicated and its variational problem 
should be treated with more delicacy. 

In Section 5, we construct approximate solutions and give energy estimates for the full system. 
Finally, we present nonlinear instability for the ideal MHD system around our three kind steady 
states in Section 6. 

We put some useful vector identities which are used throughout the paper. For any two 
vectors a and b, 

Vx (axbj =a(v-b) + (b • v) a - b (V • a) - (a ■ V) 6, (1.11) 

V (a- bj = (a- V)6+ (b ■ v) a + a x (V x bj +b x (V x a), (1.12) 

2. General variational framework 

We consider the following steady magnetic fields: 
(Case B _L g) 

B = (0,0, B (or)), 

(Case B || g) 

Bo = (Bo, 0,0). 

In the case of B \\ g, we assume v 3 = B 3 = 0. By integration by parts, using (jl.7J) . (jl.5J) . (jl.2p 
and completing the square with respect to V • v, we obtain the following decomposition for 
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< L (v) , v >: for any v = (vx, v 2 , v 3 ) , 

< L(v),v>= J J L(v)-vdxdy (2.1) 



-F (v 2x ) - G (V • v, v x , v lx ) + H (vx, v lx )] dxdy. 



D 



We state some important properties which are satisfied by the above functional: 

(1) L is variational, i.e., for any u and v 

< L(u) , v >=< u, L (v) > . 

(2) < L (v) , v > is concave with respect to vi x , v 2x , v 2 respectively. 

(3) < L (v) , v > — > — oo as v x x °o and < L (v) , v > — > — oo as v 2x — > oo if F ^ 0. 

(4) F > 0, G > 0, H (0) = when H = H (vx) alone, and F, G, H are all quadratic. 

Remark 3. G = in the incompressible case. 

Remark 4. For the case B _L g, we have F = 0, H = H (vx) , and G — G (V • v,Vx) with 
6 = while we have F = F (v 2x ) ,H = H (t> 1; Vx x ) , G — G (V • v, v x , Vi x ) in the case B \\ g as 
we can see in Section 4- 

We will show that for any fixed wave number k e N, the corresponding eigenvalue > 
for the linearized MHD system takes the variational formulation closely related to the above 
variational structure. A normal mode is of the form: 

vx{t,x,y) = vx (x) cos (ky) exp (X k t) , (2.2) 
v 2 (t, x, y) = v 2 (x) sin (ky) exp (X k t) , 
v 3 (t, x, y) = v 3 (x) cos (ky) exp (\ k t) , 

where k is a wave number. Substituting (|2.2j) into (jl.8j) yields the following second-order ODE 
for v: 

AfcPo (Vi, v 2 , v 3 ) = L (vx, v 2 , v 3 ) . (2.3) 
Note that L 3 (vx,v 2 , v 3 ) = and thus v 3 = 0. We now make the following variational formulations: 

, 2 < L(u) ,u > 

X 2 k = sup ^- , (2.4) 

U=(ui,U2) ^) 

u 1 eH K (0,2Tr) 
M 2 e-H" K ~ 1 (0,27r) 

where k = 1 or 2, d y is replaced with multiplying by k in L (u) and integrations here are with 
respect to x over (0, 27r). Indeed, 

Incompressible case: we use divergence-free condition to reduce to 

C[-F(-^)+H(u,u x )]dx 

A k — SUp 



where k = 1 or 2. 



2rr 



u2 + w 



dx 
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Compressible case: 



/o * [-F ( u 2x) - G (u lx + ku 2 , Mi, u lx ) + H («!, u lx )} dx 

SUp — ^ : 



A fc = sup 
/r«(i 

n 2 G-H' K - 1 (0,27r) 

where k = 1 or 2. 
Let 



,2_ JJ D H(v,v x )dxdy 



A J = sup ^ 7 ; y -, (2.5) 

where k = 1 or 2. Then we show A 2 is the least upper bound for {^l} keN - 
Lemma 1. 

lim \\ = A 2 . 

fc^oo 

Proof. Since F (0) = 0, it is easy to see that our Lemma is true for incompressible case by 
letting k — > oo. We now treat compressible case. 

Note that, with the choice of u 2 — — \ (aui + bui x ), we have 

2 < L(u) ,u > 
\ k = sup r 

U = (UI,U'2) \ U 7 U ) 

Mie// K (0,27r) 
u 2 eH K - 1 (0,2Tr) 

iT [~ F (~I ( aUl + bu ^)x) + H U ^)] dx 
^_ S VI p 27r 2 * 

wiei? K (o,27r) J " p [u\ + p (aui + bu lx ) ] dx 

Thus, we have the following inequality 



Jo * [~ F j-J Hi + hu i*)x) + H ( u u u ix)] dx 2 

SU P , 7^ r 2 , i / \2T~^ - A fc - 

ui6i?«(o,27r) J Q p [uf + p (awi + oui x ) J dx 



In* i~ F ( U2x ) ~ G ( Ulx + ^ M2 ' Ml > u ^)\ dx IT H Ulx ) dx 

SUp — 2^ h SUp 

Ul6 if«(o,27r) / p [w 2 + w 2 ]cfe uieff-CO,^) Jo PoK + w l]^ 

« 2 e^ K HO^vr) «2e-f/ K - 1 (0,27r) 



f 27V H (u 1: u lx )dx 



< sup 



ui6H"-i(0,27r) J Q 2W p Q uldx 

Thanks to 



— (aui + buix) — > as fc — > oo, F (0) = 0, 



letting fc — > oo yields 



J 2w H (u!,u lx ) dx 



lim A 2 , = sup 



uieif^-^o^Tr) J p u(dx 
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lim A j 



k— *oo 



sup 

u 1 eH K - 1 (D) 



JJ D H (u u u lx )dxdy 
IId Paujdxdy 



since H K is dense in H K 1 . Thus the proof is complete. | 

We now show the existence of maximizer for the variational problem (|2.4|) . Assume that 



For fixed k, let 



(U,u) 



Xl 



p2ir 

I pQ \u\ + u 2 ,] dx = 1. 
Jo 



sup 

«ieH K (0,27r) 
u 2 £H k - 1 (0,2it) 



< L(u) ,u> . 



(2-6) 



(2.7) 



Lemma 2. For any fixed k, there exists a smooth maximizer for the variational problem \2. 7| ) 
with the constraint \2.b\) . 

Proof. Let u 2 } be a maximizing sequence with the constraint (j2.6|) . Then and M2 
converge weakly in L 2 (0, 27r) to u\ and m° respectively and we have 

< L (u n ) , u n >— > A?,. (2.8) 



Case lB lg, where F = 0, H = H (m): 

Since < L(u) ,u >— > —00 as «i x — > 00 (Property 3 of the functional), u\ x is bounded in 
L 2 (0, 27r) uniformly in n. Thus there exists a weak limit {u^ 1 , u!]} such that 

^ u\ x weakly in L 2 (0, 2ix) , u 2x u° weakly in L 2 (0, 27r) 
-> w° strongly in L 2 (0, 2vr) and u\ G i/ 1 (0, 2tt) . 



■it 



Next we show that {w^, ""2) i s a maximizer and satisfies the constraint ()2.6|) . Since {w^, M2} * s a 



weak limit of {u^uV,}, (u°,u°) = Jq^Po 

Po [«f 



0N2 



+ (U 



0\2 



< 1 by lower semi-continuity of L 



2. 



u: 



;) 2 ] dx 



2n 



Po 



K) 2 + K) 



0\2 



2 ^2p [it° (v»-v°)]dx 



+ 



2tt 



0\2 



> 



/■27T 

/ 2p [«? « - u$) + u° 2 [u n 2 - ujj)] dx -> as n -> 0. 



due to the fact that {w?,^} is a weak limit of {ui,u 2 } in L 2 (0,27r). We use concavity of 
the functional < L{u) ,u > with respect to u\ x , u% x , U2 (Property 2 of the functional) and the 
strong convergence of to u\ in L 2 (0, 2tt) to deduce 

<L{u°),u°> >A 2 . (2.9) 
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Let 

J («2x, U lx , U 2 , Ml) =< L (u) , U > . 

Then by concavity of this functional J, strong convergence of u" to u\ and weak convergence 
of Ui x , u 2 , Ui to u1 x , u 2 , u\ respectively, we have 

J (u 2x , U lx , U 2 ,U 1 ) — J [U 2x , Ui x , U 2 , u l) 

= J (u 2x , u lx , u 2 ,u 1 ) — J (u 2x , u lx , u 2 , u-jj + J (u 2x , u lx , u 2 , Ui) — J (u 2x , u lx , u 2 , u ± ) 

< J «, u n lx , u n 2 , <) - J «, u n lx , u n 2 , u°) + V J «, u° lx , ul ui) ■ « - <, < - u° lx , ul - 

—>■ as n — > oo. 

Thus we obtain (J2.9)) . Notice that we can not have w° = u^} = {0, 0} a.e. since, by strong 
convergence of ii™ to u\ and by (|2.8jl . we have 

> A? > 0, 



/ H(u° 1 )dx= lim / H (u?) dx 
Jo n ^°°Jo 



where we have used the property 4 of the functional. Suppose now that (tt°,tt°) = a 2 < 1. 
Then let (m?,^) = ~ ( w i) M 2) so that (S ,^ ) = 1. By the above argument (|2.9jl . we have 

<L(u°),u°>>f >AL 

leading to a contradiction. Thus {w?,^} is a maximizer satisfying the constraint ([2.6)1 . 
Case 2 £? || g, where F (m 2x ) = -Bo^iar 

In this case, by property 3 of the functional, u 2x is also bounded in L 2 (0, 2tc) uniformly in 
n and thus both and u 2 converge strongly in L 2 (0, 2tt) to and u 2 respectively. Then we 
have 

"2tt 

dx = 1. 



Jo 1 

In a similar manner, we obtain (J2.9)) and hence u° is a maximizer. 

We finally show such a maximizer satisfies the generalized Sturm-Liouville equation (J2.3j) for 
both cases. For r G E and to = {wi, U7 2 } G H K (0, 27r) x (0, 2n) , define u (r) = u° + ™, 

then by 1)4. 8 J) , we have 

<L( M (t)), M (t)><A 2 ( W (r),u(r)). 

Set 

7 (r) =< L (ti (r)) , u (r) > -A 2 fe (« (r) , u (r)) , 
then we have J (r) < for all r G R and J (0) = 0. This implies 

J'(0) = < > + < L(m°) ,w > -2A^m°,w) 

= 2 < L (u°) - X 2 k p u°, w >= for all w, 
since L is variational (Property 1 of the functional). Thus u° satisfies a normal mode, i.e., 

\ 2 kPo u° = L (u°) . 

Since p , po, and B are smooth, w° is also smooth. This completes the proof. | 
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3. Linear growth rate A 

In this section, we show that A is the optimal growth rate for the linearized system and it 
serves as the spectral radius of the linear operator. We state a global existence of solutions to 
the linearized system, which can be obtained by a straightforward method. 

Lemma 3. There exists a global in time solution (a, v,B) £ C ([0, T] ; H s (£))) to the linearized 
MED system \TJ\) . 

Let (a, v) be a solution to (jl.7|) and let A 2 > 0, then we have 

Theorem 3. 

l|a (t) ,v(t),B (t) I < Ce At \\a (0) ,v(0),B (0) || a+2 . 



where C = C (A, p , s) and A > 0. 

Proof. We show by induction on s, the number of x-derivatives. We first treat the case s = 0: 
Multiply ()4.6|) by v t and integrate over x and y, then we have 

Jt^ Vt,Vt ^ = Tt <L ^ ,v> (3,1) 

since L is variational. Notice that using (|2.5J1 yields 

H (vx,vi x ) dxdy < A 2 (v h vi) < A 2 (v, v) . 
By integrating f|3. lj) over time and by ()2.1|) . we obtain 

(v t , v t ) + J J [F (v 2x ) + G (V • v, v 1} vi x )} dxdy 

< JJ^H (v 1} v lx ) dxdy + \\<t(0),v (0) , B (0)||J . 

< A 2 (v,v) + \\a(0),v(0),B(0)\\l, (3-2) 



and 



Thus, we have 



< \\v t \\ < \\<t(Q),v(0),B(0)\\ 1 + A\\v\\. 



lull < Ce 



At 



where C = \\a (0) , v (0) , B (0)1^ and from now on, we will just use universal constant C which 
varies, only for notational convenience. Notice that these estimates exactly apply to the t- and 
y-derivatives of any order of v and a since the variational structure of (jl.Hjl is not destroyed by 
taking t- and ^/-derivatives. 

Case 1. Bo _L g where F — 0, G — (7^0 + -B 2 ) ■ v + ^^r ^ij (G = in the incompress- 
ible case), H = H ( Vl ) , (Q • V) 5 + (S • V) Q = 0: 
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By (J3.2j) and (jl.7j) . we also have 

||V • u|| < Ce At , \\a\\ < Ce A \ \\B\\ < Ce At . 
Next, we consider the case s = 1. By taking curl w-equation of (|1.7jl . we have, with uo = V x v, 

Po (^)tt = - v Po x v tt - g<rtyk, (3.3) 

where k is the unit vector in the ^-direction. Since a ty and v tt have no x-derivatives and 
Bi x = —B 2y and so have the growth rate A as in the previous step, we have 

HI < Ce At , 

where C = \\a (0) ,v(0),B (0)|| 3 . Thanks to the identity A( = -V x (V x () + V(V • () for 
any ( = {(1,(2)1 we conclude that all the first derivatives of v have the same growth rate A. 
Now for Vcr, we use the vector identity (jl.llj) to get 

B t = V {v x B ) = -B Qx v x k - kB V ■ v. (3.4) 

Plugging (HOI) and (|T3]l into dEHJ) yields 

p «« = V (B B 0xVl ) + V ( [ 7 Po + Bl] V • «) + <r t £ (3.5) 

By induction hypotheses and cr- and -B-equations in (jl.7j) . we deduce 

II V (V • v)|| < Ce At , ||Va|| < Ce A *, ||VB|| < Ce At . 

For higher derivatives when s > 2, we use induction hypotheses. Suppose all the derivatives of 
order < s have the growth rate A. Let a be multi index whose order of x-derivative is s — 1. 
Then, by taking d a of the curl of (|4.6p . curl(9 Q f) satisfies 

Po = (V x (p u)) tt + Po (9 Q w) tt - gd a a ty k. (3.6) 

The right hand side of (J3.6|) contains derivatives of v and a whose x-order < s. Clearly the 
x-order of d a a ty is still s — 1. The first and the second term of the RHS together are of the 
form 

£ C^V x [(dppo) (d,v tt )} . (3.7) 
l^l>i 

/3+7=a 

Notice that only purely x-derivatives dpp with |/3| > 1 in (J3.7|) remain nonzero since p (x) 
depends only on x. Then (J3.6J) yields 

||V x 9 a v|| < Ce A *. 

For V ■ d a v , we take dp of the first component of (|3.5j) with dpd x = d a to get 

dp (p v lt t) = d a (BqBqxV-l) + d a ( [yp + B%] V • v) + 9^«r t ^. 

Using the induction hypotheses as for the case s = 1, we deduce 

IML- < Ce At , \\a\\- s < Ce A \ \\B\\- < Ce A \ 

where C = \\a (0) , v (0) , B (0)||- +2 . In the incompressible case, we use ()3.3|) together with the 
V • v — to deduce the Lemma. 
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Case 2. B \\ g and incompressible case where F = BqV% x , G = 0, H = H (vi, vi x ): 

From ()3.2|) and div v = and using the induction hypotheses, we obtain 

IK|| <Ce A MK|| <Ce M . 

In this case, by (ll.llj) . we have Q = B v x and B ■ Q = BqV 1x . Then (II. 8|) yields 

p v tt = B 2 v xx - V (B -Q+p t ) + a t g. (3.8) 

Taking the curl of ()3.8|) yields 

p Q u tt = B%uj xx - Vp x v tt - ga ty k. 

Using the induction hypotheses and div v = 0, we can deduce the Lemma. Therefore the proof 
is complete. | 

4. Examples 

In this section, we present examples of steady states which give rise to different results. 
(1) Magnetic field is transverse to the gravity B _L g for incompressible fluids: 

B =(0,0, £ 0*0). (4.1) 

Linear operator L reads 

L (v) := -V (p t + B -Q)- gp 0x v u 

Sturm-Liouville equation is 

A? 



Variational formulation for is 



J Q 27T L{u) ■ u dx 



K= sup o J » r v >0, (4.2) 



ft 

where 



dx 



L{u)-u =H(u) = -gp 0x u 2 
Spectral radius A > is given by 

.o ffr>H (v) dxdy 

A 2 = sup JJd 1 { y > 0. (4.3) 
^eL 2 (D) JJ D Po v dxdy 

The sufficient condition for instability in the case (1) is given by 

Pox ( x o) < for some xq G (0, 2tt) , (4.4) 
(2) Magnetic field is transverse to the gravity B _L g for compressible fluids: 

B = (0,0, B (x)). (4.5) 
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Linear operator L, Sturm-Liouville equation, and £1X6 clS follows: 

L{v) = - V (-B • Q) - V (poxVi + 7p V • v) - g (p 0x Vx + p V • v) 
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(4.6) 



**A) («i > ^2, f 3) = L (vi, v 2, u 3 ) , (4.7) 
where L = (L 1 ,L 2 ,L 3 )is a linear functional in v = (vi,v 2 ,Vs) such that 

L 1 = (gPoVi+(lPo + B%)(vi x + kv 2 )) x -g(PoxVi+Po(vix + kv 2 )), 
L 2 = -k (gp vi + (7^0 + Bq) (v 1x + kv 2 )) , 
L 3 = 0, 
and so v 3 = 0. 

, 9 fn * Liu) -u dx 

x l= su p $ rri 21 3 >0 ' ( 4 - 8 ) 

MlG-ff^O^Tr), 
u 2 Gi 2 (0,27r) 



where 



L{u)-u = -G (ui,ui x ,u 2 ) + H(u\) 

7Po + Bl 

— nrir. I 7/.? 

,7Po + #o 



G (ui, u lx , u 2 ) = - (7P0 + Bl) \ u lx + ku 2 + - — ^ o2 ni 



The spectral radius A is given by 

„ ff n H (v) dxdy 

A = SU P ff 2^/ >0 - ( 4 - 9 ) 
feL2(D) JJ D PoV z dxdy 

The sufficient condition for nonlinear instability is given by, for some x$ G (0, 2tt) , 

!j2/> " M > nPar M . (4.10) 



7Po + Bl 

Notice that density inversion is not required for instability unlike the incompressible case. Thus, 
even if light fluid is on top of heavy fluid and magnetic field lines stabilizes, such steady state 
flows can't sustain themselves against small initial disturbances under the effect of compress- 
ibility. It means that the destabilizing effect of compressibility dominates over the stabilizing 
one of magnetic field lines. Furthermore, letting 7 —>■ 00 yields exactly f)4.4j) and the instability 
criterion for the incompressible case can be recovered from compressible ones as the limiting 
case. 

(3) Magnetic field is parallel to the gravity B \\ g for incompressible fluids: 

B = (B , 0,0). (4.11) 
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In this case we consider density profile p whose gradient has negative average over (0, 2tt) , 
which is stronger than (J4.4j) : 

-gpox ( x ) dx > °) 



Linear operator L, Sturm-Liouville equation, and A& are given by: 

L (v) := BqV xx - V (pt + Bq-Q)- gp 0x v t , 
A 2 B 2 

p X 2 k vi = B 2 d lxx + -f (p vi x ) x - -j^v lxxxx - gp 0x v u (4.12) 



J Q 27r L{u) ■ u dx 



\\ = sup 



u + # 



dx 



> 0, (4.13) 



where 



L(u) ■ u = -F [-^)+H(u,Uv 



1 WXJ 1 



U ™\ ~2 U xx tt(„, „, \ _ „„ „,2 d2„,2 



/. j B oj^i H ^ u *) = -9Pox u ~ K U V 
A has the following formula: 

„ ffr,H (v, v x ) dxdy 

A 2 = sup JJd \ x) >0. (4.14) 
„ 6 jfi(r>) JJ D p v z dxdy 

Lemma 4. A 2 m d^.i^| ) positive. 

Proof. We can take a family of test periodic functions in H l which guarantees the positivity of 
A. Let c = minj, p (x) and 

x 

<f) n (x, y) = c for < x < 2tt. 

n 

Since 



BlJJ^lfdxdy 



An 2 B 2 

► as n — > oo, 

rr 



we choose large n to get positivity of the numerator of A . | 

(4) Magnetic field is parallel to the gravity Bq \\ g for compressible fluids: In this case, we 
obtain linear stability instead of instability and we have the following: 

L(v) = B\v xx - V (-B • Q) - V {p 0x vi + 7p V • v) - g(p Qx V! + p V • v) , 
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where 



Q L w L(u)-udx 
A' = sup -§ — < 0, 

u=(ui,u 2 ) L Po \ U 1 + U 2 ] dx 
u 2 £H 1 (0,2tt) 



L (u) ■ u 

F (u 2x ) 

G (ui,u lx ,u 2 ) 
H {ut,u lx ) 



-F (u 2x ) - G (ui, u lx , u 2 ) + H (in, u lx ) 



B 2 v 2 



(iPo + Bl) 
1 



2 Ui x + — — —f^ + ku 2 



JPo + B 2 7p + B 2 



Wo + B 2 



Bl 



[gp ui - B 2 u lx ] - gp 0x ul - Bl 
(gpoUt + ~?p u lx ) 2 < 0. 



2„.2 



a 



lx 



and 



A 2 



(7Po + Bl) 7p 

ff D H (v,v x ) dxdy 



veH?(D) Hd Po v2 dxd V 



< 0. 



(4.15) 



5. Energy estimate and approximate solution 



In this section, we construct approximate solutions using a method originated by Grenier in 
[Hj and we do energy estimates for the fully nonlinear MHD system in both incompressible and 
compressible cases. 

We first construct approximate solutions. In our construction, 5 > is an arbitrary small 
parameter, and 9 is a small but fixed positive constant (independent of 5). We fix k with 
A = Afc (dominant eigenvalue) so that 

A < 2A. 

We define T 5 by 

9 = 5exp(XT 5 ), (5.1) 

or equivalently, 

I*-ih,§. (5.2) 

We may write the full system (jl.6j) in vector form for w = (a (t, x, y) , v (t, x,y) ,p (t, x,y) ,B (t, x,y)): 

w + A 1 (w) d x w + A 2 (w) d y w + L(w) = F (w) . 

An approximate solution w a (t, x, y) = (cr a (t, x, y) , v a (t, x, y) , p a (t, x, y) , B a (t, x, y j) is of the 
form 

N 

(t,x,y) = ^25 J r j (t,x } y) , (5.3) 



w 
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such that 

w a + A 1 (w a ) d x w a + A 2 (w a ) d y w a + L (w a ) = F (w a ) + R a N . (5.4) 

We show existence of such approximate solutions in the following lemma. The key point is that 
we can choose a dominant eigenvalue A with A < 2A to make this construction work. 

Lemma 5. For any fixed integer N > 0, there exists an approximate solution \5. fy) satisfying 
\5.4\j - Furthermore, for every integer s > 0, there is small 6 > such that ifO<t<T s as in 
A5.ty) . the j-th coefficient rj and the remainder satisfy 

INIjj. <C(s,N) exp(j At), forl<j<N, (5.5) 
WKWh* < C (s, N) S N+1 exp {(N + 1) At} . (5.6) 

Proof. The construction of rj will be made by induction on j. The idea is as follows. We split 
the system into linear and nonlinear part: 

d t w a + A 1 (0) d x w a + A 2 (0) d y w a + L {w a ) (5.7) 

= [A 1 {w a ) - A 1 (0)] d x w a + [A 2 (u a ) - A 2 (0)] d y w a - F {w a ) =: h {5) 

with the Taylor expansion of w a in S 



N 
3=1 



Thus, our rj is the solution of the part of ()5.7|) which corresponds to the coefficient of S J in its 
Taylor expansion. 

For j = 1, take for r\ the smooth normal growing mode to the linearized system with our 
chosen wave number /c and the corresponding dominant eigenvalue A = Afc as in ()2.2|) . Clearly 
this growing mode fulfills ()5.5|) . 

Assuming that we have constructed Tj (j < N) which satisfies f!5.5|) . we construct Let 



^5 k r k (t,x,y). 

k=l 

Define the nonlinear part of the system substituted by Uj clS 

h j+1 (5) = [A 1 ( Uj ) - A 1 (0)] O x Uj + [A 2 ( Uj ) - A 2 (0)] d yUj - F ( Uj ) . 

Since this is the nonlinear part of the system and the terms in come from the terms S k for 
1 < k < j, it is enough to consider Uj in order to collect the (j + l)-th coefficient of nonlinear 
part of the expansion. Then r\,- +1 is defined to be a solution of 

d t r j+1 + A 1 (0) d x r J+l + A 2 (0) d y r j+1 + L (r j+1 ) 

_ -h™ (0) 
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with initial data r J+ i (0, x, y) = 0. Notice that 



jffi? (Q) 

O' + l)! 



jimi+j2ni2+--+jp+im p +i=i+l 



tdM m 1 m 2 r m Pf) r m P+l 



where > 0, 1 < k < j, and Bj\ depends on A 1 and F. Induction hypothesis ()5.5|) for 
rfc, 1 < k < j applies to get, for all s, 



hffi (0) 



(j + 1)! 



< C (s, N) exp [(jimi + j 2 m 2 + •• + j p+ im p+ i) At] = C (s, AT) exp [(j + 1) At] . 



Thanks to our linear estimates for A in Lemma 4 and Duhamel principle, we have 



< c 



Mt-r) 



h 



(3+1) 
3+1 



(0) 



J) 



(i + i)! 

<C(s,N) / e A( *- r) e (j+1)Ar rfr 



< C (s, N) e u+1)xt 



since j + 1 > 2 and A < 2A. 



We now define w a = Ylj=i an d ^ satisfies 



iV 



d t w a + A 1 (0) fi^ + A 2 d y w a + L (w a ) = -J2 

3=1 



(0) 



Then is defined to be the sum of all higher terms than A" in nonlinear part of the ^-expansion 



N 



R% = h(5)-J2 

3=1 



^Tif (0) 



which clearly satisfies ()5.6|) and our proof is complete. | 

We state local in time existence for the incompressible ideal MHD equations: 

Lemma 6. (Local existence to the full system) For all s > 3 and for any given initial data 
(uo, v ) G H s (D) such that p (0) = p (x) + a (x, y) > m > 0, there is a T > such that there 
exists a unique solution (ex, v,p, B) e C ([0, T] : if s (£))) to M.b]) with p (t) = p (x)+a (t, a;, y) > 
0. 

We first treat the incompressible case and we use some vector identities to get: 

a t + v ■ V (p + a) = 0, 

(p + o-) (v t + v ■ Vi>) = -V (p + + 5 ■ + B ■ VB 

+ B -VB + B- WB + ag, 



IS 



HYUNG JU HWANG 



B t = B ■ X7v - v ■ VB + B-Vv-v- V-B, 

where k is the unit vector in the z direction. Let w (t, x, y) = (a (t, x, y) , v (t, x,y) ,p (t, x,y) ,B (t, x,y)) E 
C ([0, T] ; H s (D)) be a local solution as constructed above. 

Let w a (t, x, y) = (a a (t, x, y) , v a (t, x, y) , p a (t, x, y) , B a (t, x, y j) be an approximate solution. 
We now estimate the difference 

w d (t, x,y) = w (t, x, y) - w a (t, x, y) . 

af + Va d -v + V K + p ) • v d = -R a N>1 , (5.8) 
(p + a) (v? + v d -Vv + v a - Vv d ) + a d « + v a ■ Vv a ) 

= -V (p d + ^B d (B + B a ) + B ■ B d ^j + B d -VB + B a ■ VB d (5.9) 
+ Bo ■ VB d + B d ■ VB + a d g-R a N2 , 



B d = Bo ■ Vv d - v d ■ VBo + B d ■ Vv + B a ■ Vv d - v d ■ VB - v a ■ V ' B d - R a N3 . (5.10) 

Lemma 7. For s > 3 and assume Hcr^ < ^ || Pol loo' then there exists a continuous positive 
function C depending only on s,p such that 



±\\ w d \\ 2 <C(\\w d \\ ,|K|L +1 ) \\w d \\ 2 

sjf W Us — M I lis' 11 11 s+U M \\s 

2 II T-\n I I 2 || Tin 1 1 2 



||-^iV,l|| s 11^^,2 IL II^V,3|[ S 

Proof. This energy estimate is straightforward and thus we give a brief sketch. We take d a 
of the equations (|5.8|) - (|5.1U|) . multiply through ^-derivatives, and integrate over D. For sym- 
metric terms such as Vd a o~ d ■ v in the first equation vanish upon integration due to the di- 
vergence free condition for v , v a , v d , B, B a , B d . We now estimate nonsymmetric terms. Terms 
which need attention are (B d ■ V) d a B d , (B a ■ V) d a B d , (B ■ V) d a B d from t>-equation ()5.9|) and 
[B d ■ V) d a v d , (B a ■ V) d a v d , (Bo ■ V) d a B d from 5-equation ()5.10|) . Upon integration, the three 
corresponding pairs exactly cancel out. To see this, 

B d d t d a B d d a v d = diBfd a B*d a v* - Bfdid a v*d a B* = -B d d t d a v d d a B d , 

Btd t d a B d d a v d = diB?d a Bfd a v$ - B?did a vfd a Bf = -BfdfavfcBf, 

B Ad a B d d a v d = diBoidaBp^ - Bo l d l d a v d d a B d = -Bo l d l d a v d d a B d , 

thanks to divergence free condition for v and B. Therefore we obtain the Lemma. | 

We now extend the energy estimates to the compressible case using symmetrizer, which is 
necessary in this compressible case. We write the full system (jl.6|) in components: 

o,t + ViCTj + (p + a) v iti + Po x v x = 0, 

(Po + a ) i v j,t + ViV j:i } + (p + a) qa :j + B^B^ - BiB j:i + B B 3;j 
= [PoQ (Po) - (Po + a ) l] Po,j + B i B o,i - B 0lj B 3 
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B jjt - BiV j:i + v^iBj + ViB j:i + v i:i B = -ViB 0ji . 
Here j = 1,2,3 and twice i means the sum over i — 1,2, 3. 

We rewrite the full system near the steady state ^p , 0, B j in vector 
(<j,v 1 ,v 2 ,v 3 ,B 1 ,B 2 ,B 3 ): 

w t + A 1 (w) d x w + A 2 (w) d y w + L(w) = F (w) , 
Introducing and multiplying the symmetrizer with q (p) = 7p 7 ~ 2 

D = diag (q (p + a) , p + a, p + a, p + a, 1, 1, 1) , 

leads to the following symmetric matrices 



DAi 



notations for w 



(5.11) 
(5.12) 



DAo 



DA, 



fl ( Or, -4- (T 1 7)1 


(Po + a ) Q (Po + ff ) 





o 


o 





( Dr, -4- (T I (7 I /}r> -I- (T I 

IM) ~ u } H KrO ~ u ) 


(Po + CT ) 


ui 





o 


Bo 


B 3 + Bo 








(Po + ff ) «1 





- 


-B 1 














(Po + ff ) «1 








-B l 














Vi 











B 2 


-B 1 








Vi 








B 3 + B 


-B 1 








^i 


q (Pq + o) v 2 





(Po + 0") 9 (Po + ff ) 

















(P + ff ) ^2 








-B 2 








(Po + °) 1 (Po + a ) 





(P + ff ) ^2 





B l 





B 3 + S 











(P + O) ^2 








-s 2 





-B 2 


B 1 





^2 























v 2 











-£>3 + Bq 


-B 2 








v 2 


Q (Po + ff ) V 3 





(Po + 


o) q (po + a ) 








o ■ 







(Po + 0") ^3 








-B 3 
















(PO + ff ) ^3 








-B 


3 




(Po + °) Q (Po + °") 





(p + a) v 3 


#1 


B 2 


Bo 


j 





-B 3 





Si 


^3 























^3 
















B 








^3 . 





DL (u) 



PoxQ (Po) vi 

Q (Po) P0x° + PO?' (Po) P0x° + Ate #3 



— B 0x Bi 
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PoxQ (Po) «1 - PoxQ (Po + °) Vl 

[(o- + Po) 9 (o - + Po) - Po<? (Po) - 9 (Po) Pox^ - Po?' (Po) P0x°] Pox 



/)/•'(//') 





Now that DA 1 , DA 2 and DA 3 are symmetric, we have the following local in time solution to 
the full system via standard energy estimates for symmtrizable hyperbolic system as in |3| , ^3] 
in the absence of magnetic fields. 

Lemma 8. For all s > 3 and for any given initial data (&o,vq) G H s (D) such that p (0) = 
Po (x) + o"o (x, y) > m > 0, there is a T > such that there exists a unique solution u = 
(a, v,B) EC ([0, T] : H s (£>)) to fOP wift p (t) = p O) + °" (*> ^ v) > °- 

Since construction of approximate solutions is similar to the incompressible case as in Lemma 
5 for a hyperbolic system, we omit it. We now estimate the difference of an exact solution and 
an approximate solution. Let w (t, x, y) = (er (t, x, y) , v (t, x,y) ,B (t, x, y)) G C ([0, T] : H s (£))) 
be an exact solution and w a (t, x, y) = (o~ a (t, x, y) , v a (t, x, y) , B a (t, x, y)) be an approximate 
solution as constructed in Lemma 5. Then their difference is 



w 



w — w 



[a 



a a ,v 



v a ,B-B a ), 



and it satisfies 



A* (w a + w d ) d iW d + [ Al i w<1 + wd ) ~ Al (™°)] d ^ a 



(5.13) 



i=i 



i=l 

= -L(w d ) +F(w)-F(w a )-R a N . 

This symmetrizable hyperbolic system for w d allows the following energy estimates. The proof 
is straightforward by classical energy methods as in j^j, [S]. 

Lemma 9. For any s > 3, there exists a continuous function g s ( 



such that 







|^|||'<^(|||^||| s ,llkilL + 0|||^|||' + iiii?^ 1112 



h N\\\s J 



dt 

where \ \\-\\\ s is defined by 

IIMlis = daVD (po + °") daV - 

\ct\<s 

Here D is the symmetrizer as in (|5.12j) and notice that |||-||| is related to the usual norm 

INI*. b y 

v \H H s < \\\v\\\ s < c s {\\<t\\ h . , HpqIU) \\v\\ h . , 

since p > c > and so D > rjld for some rj and s > 3. 

Notice that all three norms ||-|| s , |||-||| s , and \\-\\ H s are equivalent since p is smooth with a 
positive minimum. 
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6. Nonlinear instability 

Proof of Theorem 2. Let w a (t, x, y) = (a a (t, x, y) , v a (t, x, y) , p a (t, x, y) , B a (t, x, y)) be an ap- 
proximate solution with N to be determined later. For any 5 > 0, there exists a local in time 
solution w 5 (t, x, y) = (a s (t, x, y) , v 5 (t, x, y) , p s (t, x, y) , B 5 (t, x, y)) with initial data w a (0) to 
the full system By the Lemma 7 and Lemma 9, we have 

(/ ' ! d]]2 <C {\\w d \\ s , \\w a \\ s+l ) \\w d \\] + C^+V^^*, (6.1) 



— \\w 
dt 



with w d (0) = 0. 
Let 



T = sup{t| \\w a \\ s+1 <Ld/2, \\w d \\ g < uo/2] 



where uj is a small positive number which assures local existence. Then T depending on 5 is 
well defined since w d (0) = and \\w a (0)|| s+1 = O (5). We claim that the instability time T s 
occurs within the existence time T, that is, T s < T. Suppose not, i.e., T < T 5 . Then for t <T, 
by our construction of approximate solution, we have 

N N N 

i=i i=i i=i 

if 6* is small. Now we appeal to the definition of T and (JbM)) to get, for t < T, 

| IKH* <C(p/2,u/2) \\w d \\ 2 s + C5 2 ^e 2 ^ xt . 
We choose iV > satisfying 

C(u/2,oj/2) <2(N + 1), 
so that has growth rate at most (N + 1) A. Then using Gronwall inequality leads to 

\\w d \\ s < C5 2(N+1) e 2 ^ xt = C9 2(N+ V < u/2, 
if uj is small. Thus we can deduce T s < T. Now at the instability time T 5 , 

N 



\w a {T 5 )\\ L1 >5\\ ri {T S )\\ L1 -J2 Sj h( T5 ) 



L 1 



i=2 

N 



>C5e XTS -J^C^e^ 



J'=2 

N 



3=1 
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We then deduce at time T s , 

IKML>KML-IIK 

> \\w a (T 5 )\\ L1 - \\(w s 

>^e- ce N+1 > -e 

- 2 ~ 4 

I 
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